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The spin structure of the system of quasifree fermions having total angular momentum J = 1/2 
is studied in a consistently covariant approach. Within this model the relations between the spin 
functions are obtained. Their particular cases are the sum rules Wanzura - Wilczek, Efremov - 
(N ■ Leader - Teryaev, Burkhardt - Cottingham and also the expression for the Wanzura - Wilczek twist 

2 term gY^ ■ With the use of the proton valence quark distributions as an input, the corresponding 
spin functions are obtained. The resulting structure functions gi and gi are well compatible with 
the experimental data. 
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CNJ , The nucleoli structure functions, both unpolarizcd and polarized, are the basic tools for understanding the nucleon 
internal structure in the language of the QCD. Precision measurements on the polarized structure functions have been 
completed only in the recent years [jjj - ||. These functions contain information which is a crucial complement to 
the structure functions obtained in the unpolarized deep inelastic scattering (DIS) experiments. At the same time 
the interpretation and understanding of polarized structure functions seem to be much more difficult than in the case 
of unpolarized ones. Actually, until now it is not clear how the nucleon spin is generated from the spins and orbital 
momenta of the quarks and gluons. For the present status and perspectives of the nucleon spin physics see M and 
t-H , citations therein. The more formal aspects of the polarized DIS are explained in [jl0| and JTT[ . 

The spin in general is a very delicate quantity, which requires correspondingly precise treatment. It has been argued, 
that for correct evaluation of the quark contribution to the nucleon spin it is necessary to take properly into account 



o 

r"| ■ the intrinsic quark motion |l2| - 22 1. The necessity of the covariant formulation of the quark - parton model (QPM) 
CIh! for the spin functions has been pointed out in p3| ]. These requirements are not satisfied in the standard formulation 
i of the QPM, which is currently used for analysis and interpretation of the experimental data. 

| In the paper p2| ] we have demonstrated the role of the intrinsic motion for the spin structure functions, using a 
simple model of the system quasifree fermions on mass shell. The basic requirement was a consistently covariant 
formulation of the task for the system of fermions, which are not static, being characterized by some momenta 
• i-H , distribution in the frame of their center of mass. In the present paper we attempt to further develop this approach. 
In section || we are introducing the spin structure functions g\(x) 1 g2{x), the longitudinal and transversal net spin 
density distributions sr,(x), st(x), and the density of total angular momentum j(x). Then in Sec. |l| we show how 



all of these functions are mutually related. Finally in Sec. IV we apply the suggested approach, with some simplified 



assumptions, to the description of the proton spin structure and make a comparison with the experimental data on 
the gi(x) and (72 (x). The last section is devoted to the short summary and conclusion. 



II. SPIN STRUCTURE FUNCTIONS AND SPIN DISTRIBUTIONS 



In the previous paper [£2j we have shown that the spin structure functions, related to the spherically symmetric 
target consisting of the three quasifree fermions (of spin 1/2), having resulting total angular momentum J = 1/2, can 
be written as: 
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v 

|q| 



pi 



\q\ 2 P0 



1 - 



p|/2 
Po + m 



pois + pi |q| 
Mv 



d 3 p 
Po ' 



(1) 



*to be published in Phys. Rev. D 



1 
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where po = v m2 + P 2 j Pt = P2 +-P3> an d i s the charge weighted distribution 

3 

/c=i 

This distribution is constructed from the polarized distributions of individual fermions 

AG fc (po) = Gfc,+i/2(po) - G fe ,_i/ 2 (po), (4) 

which satisfy 

Gk{po)d 3 p = 1; Gfe(po) = Gfc !+ i/ 2 (po) + G fc _i/ 2 (po), (5) 

3 

/ AG(p )d 3 P = 1; AG(po) = J! AG fc (po). (6) 
J k=i 

Distributions Gk,x{po) measure probability to find a fermion in the state 

u 0, An) = ( P f A " ] ; -ncr0 An = A0 An , A = ±i (7) 



where the direction n coincides with the direction of target polarization J and a standard normalization is used: 

N = ^^, 4^ = 1. (8) 
Po + m 

Now, let us try using the same system to calculate some spin distribution functions. In the first step we shall 
define these distributions in terms of the fermion momenta related to the target rest frame, then we shall show their 
representation in the variable x. 

The net spin density corresponding to the projection on the direction n' is defined as 

5(p,n,n')=^G A (p )« t b,An)(n'S) M (p,An); E = | (S ff) ' (9) 
One can verify, that this expression can be modified 

t pa pn' 



S(p, n, n') = ^ E G ^°) ( Ann ' - Ann '— - 
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n2 "T" ^An , 

to) (po + m) 



(10) 



We assume, that the beam direction is defined by the vector k = (|k| , 0, 0), then one can obtain the following particular 
cases of the distribution (||). 

1) Longitudinal polarization in longitudinally polarized target, i.e. n = n' = (1,0,0), then 

3 

per • pn' = ^2p 2 i a i n' i + '^2p i p 1 a i n' :j = p\a\ + + PzPidz (11) 



and the relation ( |l0| ) can be simplified 



Sl (p) = ^- AG(po) fro + — ^— ) . (12) 



2p V Pa + m 

2) Transverse polarization in transversely polarized target, i.e. n = n' = (0, 1,0), then 
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per pn = '^pfaiTii + ^PiP j cr i n : j = p\o 2 + + (13) 



and the relation (llOj) can be simplified 



Sr(p) = 7^AG(p ) fm+ — ^— ) . i 1.4) 



2po V Po+m / 

3) In a similar way, one can also obtain the polarizations S'J(p) and S^(p), which are related to the density of 
longitudinal polarization in the transversely polarized target and vice versa. The density Sj(p) can be obtained from 
the relation ( |ic| ) after inserting n =(0,1,0), n'=(l,0,0) and S^fp) with n =(1,0,0), n'=(0, 1,0) correspondingly. 
After some calculation similar to that for obtaining the relations ([12]) and (Q), one gets 

Sl(p) = S£(p) = -LaG(po)^^. (15) 
2po Po+m 



J(p,n,n') = ^G A (po)u t (p,An) (n'j) w(p, An); j fe = S A + Z fe = ~ ( ^ & ^ ' ) - ie fcim p/ . (16) 



4) The density of total angular momentum can be defined as 

2 V a k J dp m 

One can verify, that after some calculation this expression can be simplified: 

J(p,n ) n') = ^nn'AG(p ). (17) 

This result implies, that J has rotational symmetry, so there is no distinction between longitudinal and transversal 
density: 

Jl(p) = Mp) = ^AG(po) = J(p), Jj(p) = J$(p) = 0. (18) 

Further, is it possible to express the obtained distributions as the functions of x instead of pi For a simplification 
we shall from now on assume that 

Q 2 > AM 2 x 2 , (19) 

which implies 

M = v /l + 4M 2 a; 2 /Q 2 -» 1. (20) 
Then the <5 function term, which defines the transformation p — ► x in Eqs. (g), (|) will be simplified 

in this limit the coordinate p\ defines the beam direction. The spin structure functions ([TJ), (|^) are now simplified 
accordingly: 



9x{x) = i J H(po) (m + px 
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Apparently, the convolution denned by the S function (^l|) also gives the rule for transformation of the spin distributions 
expressed in the variable p to the corresponding representation in the variable x: 



j(x) = / J(P)S -x)d 3 P=lj AG(P0)« - *) (25) 



•*(*) = / ( ^ - A *P = i I AG(,o) (m + -A-) 8 (26) 



M J 2 J ' \ po + mj \ M J p 

m.) - / - .) A- - J / (» + ^_) * (*±* - .) £. 07, 

In the last integral we could replace p\ by p^/2 because of the axial symmetry. Further, it is obvious, that the 
densities Sj and S^, expressed in the variable x, vanish due to the symmetry: 

s L (x) = s T (x) = - / AG po — ■—& — — x ) =0. (28) 

2 J po + m \ M J po 



What is the meaning of the integrals in the relations (|22j) - fl27[)? To simplify this question, let us assume the 
same shape of the distributions Gk(pa) for all the three fermions. Then the distributions AG and H differ only by a 
constant factor, in which charges and polarizations of individual fermions are absorbed: 

3 

H( Po ) = kAG(po); K = Y, e * A G k (p )d 3 p. (29) 

k=l J 

Now, in agreement with the results obtained in p2] ], one can observe the following. The relation ( p2| ) can be 
rewritten 

9\{x) = | J AG(po) (m+pi 
and after integration over x one gets 



Pl ^ S (l*±J*-l\i* (30) 



po+mj \ M J po 



Fx = J 9x(x)dx = lj AG(po) Q + ^) d 3 p = k (nE) , (31) 

where n is the direction of the target polarization and (nS) represents the resulting projection of the spins coming 
from the individual fermions. In a similar way one gets: 

j(x)dx = i J AG( Po )d 3 p = i (32) 

i.e. the integral represents the resulting projection of the total angular momentum. Further, one can easily check: 

J s L (x)dx = J s T (x)dx = I J AG(p ) Q + d 3 P = ^. (33) 

Moreover, the following relation is valid: 

K- s T (x) = gi(x) + g 2 (x). (34) 

Let us note, despite our assumption that the target consists of just the three fermions, the suggested approach 
is more general. Since the spin functions are always based on the differences like (|]), all the resulting relations are 
equally valid for any target consisting of the fermions with AGfc(po) ^ 0; k = 1, 2, 3, which are embedded in another 
system with compensated spins: AGk{po) = 0; k = 4, 5, 6, ... 

Now we can summarize: 

i) The function j(x) measures the contribution of total angular momenta (spin + orbital momentum) of the 
constituent fermions to the target spin. 
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ii) The functions S/,m(o;) measure the net spin contribution of the constituent fermions to the spin of target with 
longitudinal (transversal) polarization. Apparently, one can calculate also the corresponding densities of the orbital 
momentum as 



h(T) =j(x) - S L{T )(x). 



(35) 



Hi) Obviously the functions g%(x), j(x), sl(x) and St{x) are equivalent in the case of the static fermions, where 
orbital momentum does not play any role. The considered distribution functions AG and H have rotational symmetry 
in the target rest frame, which is a necessary condition for the target with spin J = 1/2. It follows, that meaning 
of our function j(x) suggested above does not depend on the orientation of the target polarization (longitudinal or 
transversal) with respect to the beam direction. Let us point out, the last statement can be deduced only in the 
framework of the rclativistically covariant description, in which the rotational symmetry of the target is properly 
taken into account. At the same time, let us note that in general, the function gi(x) is not equivalent to the measure 
of the longitudinal spin density sl (x) , only their integrals over x are equal (up to factor k) . 



III. RELATIONSHIP AMONG THE SPIN FUNCTIONS 



Before next discussion we shall first prove the following proposition: 
The functions V n (x) defined as 



V„(x) 



satisfy 
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(36) 



(37) 



for any powers j, k and function H, for which the integral (|36|) exists. Proof of the last relation is given in the 
Appendix A. 



A. Spin structure functions g\(x),g2(x) 



With the use of the relations (|36[),(|37|), as shown in the Appendix B, one can rewrite the Eqs. (E3J) and (124) as 



92{x) 



a(x)V (x)+ b{x,y)V {v)dy 



(38) 



where 



gi(x)+g 2 (x) = - 



A(x)V (x)+ / B(x,y)V (y)dy 



(39) 
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B(x,y) 



(x + X ) ^ X 2 — Xq 
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9 9 
XQ 



Now, one can easily check, that in the limit xq —* the relations (|38|) , (39) are simplified: 

g 2 (x) = -V (x) + f U— A - V (y)dy, 

j x v y y j 



(41) 



(42) 
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These relations imply 



and 



gi(x)+g 2 (x) 



gi(x) = V (x) 



[9i{%) +92(x)]' 

The combining of the last two relations gives 



1 / r 2 r 



A X ^ \ 1/ / \ j V {x) 

4-t - — Vo{y)dy . 

y y J x 



[ gi ( X )+g 2 ( X )]' = - 9 -M 



(43) 



(44) 



(45) 



(46) 



g%{x) = -gi(x) + j ^^-dy, 



(47) 



which is the known expression for Wanzura - Wilczek twist-2 term for 172 approximation p4| . 

Can we now obtain a similar relation for the case xq > 0, i.e. for the massive fermions? Let us combine Eqs. (|3"§|). 
(pit) to the form 



9l[x) + 92{X) = \[ { B{X ' V) ~ ^ HX ' V) ) Vo{y)dy 



and let us try to express differentiation of r.h.s. as a combination of g\ and g 2 : 

ci{x)g 1 (x) + c 2 {x)g 2 (x). 



1 ' A{x) 

B ( x i v) — rr b ( x ' v) ) v o(y) d v 



a(x) 



(48) 



(49) 



In the Appendix C it is shown, that after inserting g\, g 2 from Eqs. (|38|), ( ]3S| ) this equation is solvable for Ci(x), c 2 (x), 
then after comparing with Eq. (Est) we get 



gi(x)+ 1 + 



A(x) 



where 



ci(x) = -- 



a(x) 

X 2 + AXXQ + Xq 



92{x) 



a{x)gi{x) +c 2 (x)g 2 (x), 



(x 2 ~xl){x + 2x () y 
In the same Appendix it is shown, that Eq. (p(t) implies 



c 2 (x) 



x (x 2 + xx + 4Xq) 
(x - Xq) (x 2 - Xq) (x + 2xo) 



x-xq , . x(x + 2x ) f 1 y 2 - x% 
g 2 {x) = gi(x) + - -j- / 3 — gi{y)dy, 

X [X + Xq) Jx V 



(50) 



(51) 



(52) 



x x ~\~ 2xq r 
gi(x) = g 2 (x) 5 * I g 2 (y)dy. (53) 



X Xq X J j. 

One can check, that for xq — > both the last relations are equivalent to Eq. (fl6|). 

Further, in the limit xo — > one can also easily calculate the momenta of the spin structure functions gi,g 2 . If we 
define 

(x a ) = f x a V (x)dx, (54) 
Jo 



G 



then after integrating by parts the following relation is obtained: 



Application of this relation in Eqs. 



o Jx v p 0L + \ 

(H) then gives 



a + 1 



x a [ gi {x)+g 2 {x)]dx = {x«) 
o {a + 2) (a + 3) 



(55) 



(56) 



1 a I \ j / a\ a (a + 1) 

x g2\x)dx = — (x ) 



io (a + 2) (a + 3) 

for any a, for which the integrals exist. The last two relations imply 

1. 1 



.T 



a + 1 



Si(z) +ff2(x) 



dx = 0, 



(57) 



(58) 



which for a = 2,4, 6, ... corresponds to the Wanzura - Wilczek sum rules p4| . Other special cases correspond to the 
Burkhardt - Cottingham (a = 0) and the Efremov - Leader - Teryaev (ELT, a = 1) sum rules [p5|, pq|. 



B. Spin distributions j(x), sl(x) and st(x) 

Now, let us try to find the relations among the spin functions j (x) . Sl(x), St(x) and the structure functions g±, g 2 . 
For this purpose we only slightly change the definitions (|2^),(p6|) ( p7j ) in which we replace the function AG(po) [Eq. 
(Q), the sum of spin contributions from individual fermions] by the function H(p ) [Eq. (H|), spin contributions are 
weighted by the charges, see also the paragraph involving Eq. (po])]. The new relations read: 

j(x) = \[ H( PQ )S ( P°±?l - x) d?p, (59) 



2 / V M 



M*) = \ I H(Po) (m + J-) i (*±» - x) ^, (60) 
2 7 V po + mj \ M J p a 



s t(x) = ~J H(p ) 
Now, using a standard notation 



P 2 T \s(P°+Pl x \ fP (61) 



2(p + m)J \ M J pa 



9t(x) =gi(x) +g%{x), (62) 
we get from the relations (|24|) and ([n]) the equivalence 

g T (x) = s T (x). (63) 
Comparison of the relations (|3^) and (^9|) implies, that 

j(x) = ^Vo(x). (64) 

In the Appendix D we have shown, that the distribution Va(y) can be extracted from the relations d38|), (|39]), so we 
get: 

x 2 +^ 3x 2 + 2xx + 3x 2 f 1 y 2 -x 2 g f 1 ( x(y + x ) 2 \ y 2 - x 2 

3{x) = ~Yr*~ 9l{x) + — or j. \ 2 — / —^—9i{v)dy + / In 2 —r—gi{y)dy, (65) 

^ 2(x + x ) y Ji u(i + io) / y 
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x 2 +xl 1 f 1 6y 2 x - 2yx 2 + x 2 x + 2yx\ - 
j( x ) = -7T-7I — TT^W - 7773 37/ 3 92{y)dy. (66) 



^(x-Zor v ' 2(x 2 -^)j; y 2 
Then for xq — ► we have: 

y\ 9i(y) 



j(x) = \gx{x) +\ J (3 + 2 In: 



x/ y 



dy, (67) 



j(x) = -^(i) - I /' (31/ - x) ^Mdy. (68) 



Further, comparison of the relations (|6^) and ( f22| ) implies 

Si (x) = ft («) - i J H( P0 ) Pl S (»±£ - *) tE. (69) 

The last integral can be expressed by means of the function Vo (see Appendix E): 

1 x 2 +x 2 x-3x f 1 y 2 - xl f 1 ( x(y + x ) 2 \ y 2 - x 2 

SL[X) = 2-^ 9l(x) -WT^l —» MIU,!I - J, ln ( ) — 9l{y)dy - (70) 



A similar procedure gives also 



1 x 2 + xl f 1 / 2y-x x + 2x 

2x(x-x ) 92[X)+ J x \ 2y 2 x 2 -x 2 



= + I ( ~ ) 9Mdy (71) 



One can easily check, that for xq — > it follows 



dy, (72) 



sl(x) =-lg 2 {x)+j (~-^Jg 2 (y)dy. (73) 

Finally, combining relations (|52]),(p3|) and ( p3| ) one gets 

r \ Xq / \ x (g + 2x ) f 1 y 2 -^ , \ j (1 a \ 

g T [x) = —gi(x) + — -y- / 5 — gi{y)dy, (74) 

* (z + xo) Ji 2/ 

9t{x) = g 2 (x) - ^ ? + / g2(y)dy. (75) 

and for — > it follows 



y x 

In the end of this section let us point out, that the simple relations above, which define mutual transformations among 
the functions gi(x), g2(x), j(x), Sl{x) and st(x) are obtained on the assumption that the fermions have some fixed 
effective mass xq. In a more general case, for example when the structure functions are related to the system of 
fermions with some effective mass spectrum UM, such simple transformations do not exist. 



X Xq X *^Q J x 



m (■> )= I — </;/ - [ 1 g2(y)dy. (76) 
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IV. VALENCE QUARKS 



Now let us try to apply the suggested approach to the description of the proton spin structure. For simplicity, as 
in (||], we assume: 

1) Spin contribution from the sea of quark-antiquark pairs and gluons can be neglected. Then the three fermions 
in our approach correspond to the three proton valence quarks. So, in this scenario, the proton spin is generated only 
by the valence quarks. Let us remark, in the cited paper we have shown that this assumption does not contradict the 
experimental value T±. 

2) In accordance with the non-relativistic SU(6) approach, the spin contribution of individual valence terms is given 

as 

s u = 4/3, s d = -1/3. (77) 

Let us denote momenta distributions of the valence quarks in the target rest frame by symbols h u and with the 
normalization 



2 

then the generic distribution (m reads 



I I h u ( Po )d i P = I h d (p )d a p = 1, (78) 



1\ 2 1 



H(po)=J2 e * Ah i^= U) 3 hu ^-{3j t dijK)) - (79) 

In the papers fll9| , using a similar approach, we have studied also the unpolarized structure functions. In 
particular we have suggested, that the structure function F2, can be in the limit (119), expressed as 

F 2 (x) = x 2 f K( PQ )-S (^7^ - x) d 3 P ; K(p ) = ]T e 2 h g (p ), (80) 
J Po \ J q 

where h q are distributions of the quarks with charges e q . For the valence quarks one can write 

4 1 

Fi{x) = gXu v (x) + -xd v (x), (81) 

then the Eq. ( |S0| ) can be split: 

q v (x) = x J h q (p )^-5 (^ P °^ I Pl - rf3 P; q = u,d. (82) 

In an accordance with the definition (^), in which h q is inserted instead of H, one can write 

q v {x)=xVl l {x), (83) 

then the relation (|37]) implies 

VS{x) = \ (x + ^ ) Vl^x) + \ £ (l - ^ Vl^dy, (84) 



which after inserting from Eq. (p3[) gives 



V q (x) = \ 



l + X 4)q V (x) + f *M*y 

x J Jx V y J y 



(85) 



Obviously, the function Vo (x) generated by distribution (|79|) according to the definition ( p6|) can be decomposed 

Vo(x) = ^V u (x)-±V d (x) 

and if we define 
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w g {x) = —u v (x) - —dv(x), 



(87) 



then one can check, that inserting Vq from the relation ( |S6| ) to the relations (J3q) , ( |39| ) with the use of Eqs. (p5[), ( |87j ) 
gives: 



u> fl (a;) — 2 (a; + x ) z 



2 j/-x w g (y) 



dy 



(88) 



02 00 



II )w g {x)+3(x + x ) / -3-2 

x > J x (y + x y y 



dy 



Obviously, the structure functions can be split into the two parts, corresponding to u and d quarks 



2\ 2 2 



1\ 2 1 



9i(*) = [ 1 1 l$[x) - ( § J 0=); J = 1, 2, 



where the partial structure functions read: 



9 q x(x) 



/ x / . n2 f 1 y-j^qviy), 

q v {x) - 2{x + x ) / -3 dy 



(89) 



(90) 



(91) 



02 0) 



1 



ft x °\ / \ 1 o f 1 n2 l 1 y-^^qviy), 

\1 - — ) qv(x) + 3(x + x ) I - — 3 — - — dy 



x I 



x {y + x ) y 



q = u,d. 



(92) 



Now we can express the corresponding contributions of different quarks to the spin distribution functions. Apparently 

(93) 



40*0 =9i(x) +gl{x) = 



xo , s , , , ,2 f y-xp qy(y) 

—q v (x) + (x + x ) / — 73— —dy 

x J x {y + x ) y 



Further, after inserting from the relations (pi]), ( |8^ ) and ( |S3"| ) to the Eq. (El) one easily gets 

%o\ iv{y) 



4(*) = 7 



t 1 x o\ 1 \ a i , \2 f 1 y-xo iv(y), 

l + -t)qv{x)-4:{x + xo) / 3 dy 

x / Jx (y + xo) y 



1 _ 

2 



y J y 



dy 



(94) 



Let us remark, the Eq. (El) is obtained from the generic distribution H(po), in its place we now have the distribution 
h q (po). Similarly, a comparison of the relations ( p4[ ) and ( |85| ) gives 



fix) = - 



(95) 



Now, the net complete spin distributions can be obtained by adding individual valence terms with the weights fl77|), 
taking into account their normalization (k8|). If we define 



2 1 

w s (x) = -u v (x) - -d v (x), 



(96) 



then the complete spin distributions can be obtained from the relations ( |93| ) - (95), in which the distribution qy is 
replaced by w s . Then for xq — > we obtain: 



9i(x) = 2 



w g (x) — 2x 2 



^dy 



(97) 



92 (x) 



-Wg(x) + ix 1 



w g (y) 



dy 
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s T (x) = ^x 2 J 



i 2 f 1 w s(y) 



dy, 



(99) 



sl{x) 



w s {y) 



dy 



(100) 



3{x) 



w s {x) + I — - — dy 



y 



Let us note, if one assumes uy(x) ~ 2dy(x), then the following substitution can be used: 

/ \ 5 F2 V al , > Flval 
w 9\ x ) ~ n ' Ws \ X > ~ • 

9 x x 



(101) 



(102) 



Further, let us make a remark to the normalization of the distributions above. To simplify this consideration, we 
assume the case xq — > 0. The relation implies 



Since the relation ([57]) implies 



then using the relation (D7) one gets 



i ,i 
qv(x)dx = / xV^ 1 (x)dx. 
Jo 



Vl 1 {x)= 2 -V«{x)-2 f^dy, 
x J x y 2 



i ,i 

xVl 1 {x)dx= / V q (x)dx. 
ii Jo 



(103) 



(104) 



(105) 



From the definition 



V q (x) = h q ( P0 )S 



,1 - / /, I A / P0+Pl x ) ,/' /( 



(106) 



one obtains 



V q (x)dx = / h q ( Po )d 3 p 



(107) 



This relation combined with (105) and (103) gives 



q v (x)dx = / h q (p a )d p, 



(108) 



which in an accordance with the normalization (|78[) implies 



1 



Uv(x)dx — / dy(x)(ia; = 1. 



Now, one can also check the normalization of the functions (|97]) - ( |l0l| ). Taking into account, that 



(109) 



w s (x)dx = 1, / w g (x)dx = — , 
ii Jo 9 



then after integration with the use of relation (D7) one gets 



(110) 



j(x)dx=\ ) , I s T (x)dx = / s L (ir)efe = i 
^ Jo Jo o 



(111) 
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J^ 9l {x)dx = ^ J^g 2 (x)dx = 0. (112) 



The meaning of these integrals was discussed in Sec. ||. The first integral ([D7]) represents the sum rule on the total 
proton angular momentum J = I + s = 1/2, which has been discussed e.g. in p7[ . The integrals on the net spin 
contributions st, sl are correlated with the T\, which reaches its minimal value for xq — > 0, as we have discussed in 
p2fl , see also Eq. ( |3l| ) above. 

So, the obtained formulas enable us to calculate the spin functions from the input, in which only the valence 
distributions are used. For simplicity we shall now consider only massless quarks (xq — > 0) and for the valence 
functions xUjAx) and xdy(x) we use the parameterization obtained (for Q 2 = 4GeV 2 /c 2 ) by the standard global 
analysis in Pq| . In the Fig. [l]a we have shown the result of our calculation for g\ according to Eq. (^) together 
with the experimental data represented by the new parameterization of the world data on g\ Q for Q 2 = AGeV 2 /c 2 . 
The calculation qualitatively agrees well with the data, however it is apparent, that the data are above our curve. 
This can be connected first of all with our simplification for xq — > 0, when the Ti is minimal. In accordance with Eq. 
(|ll2|) we obtain Ti = 0.093, but experimentally rf = 0.118 ± 0.004{stat.) ± 0.007(syst.) at Q 2 = 5GeV 2 /c 2 [|. Just 



this difference is exposed in the figure. Further, in the Fig. we have shown the gi according to Eq. (|9^) and the 
precision measurement recently published by the E155 Collaboration ||. The agreement with the data is very good 
and it can suggest, that dependence of the function gi on the mass terms is rather weak. At least in our approach 
Ti does depend on mass, but T2 = regardless of the mass. In the Fig. || the corresponding spin distributions 
j, st, Si are shown for whole proton and also separately for u and d valence quarks corresponding to the assumed 
SU(6) symmetry, which gives the fractions (|77|). Figure [l]a and the left part of Fig. || also demonstrate, that the gi(x) 
and sl(x) are not equivalent. Slightly different shape of the distributions on st,sl is due to variable x, in which 
longitudinal and transversal (in respect to the beam) quark momentum components are not involved in a symmetric 
way. Otherwise, for given direction of the proton polarization, quark spin density cannot depend on the direction, in 
which the probing beam is coming. 

In the end, let us remark, that another possible effect, which can in our approach contribute to an underestimation 
of Ti is connected with the assumption (|77]). For example, if one assumes full spin alignment of the u valence quarks, 
then 

s u = 2, s d = -1 (113) 
and instead of the generic distribution ([87]) one gets 

4 1 f 1 7 

w g (x) = -u v (x) - -d v (x), J w g (x)dx = -, (114) 

which implies Ti = 7/54 = 0.13. Obviously, assuming isotopic symm etry , the same procedure can also be used for the 
neutron. For the isotopic count erpa rts of the compositions ([77]) and (113) one gets T™ = and T™ = —1/27 = —0.037, 



respectively. The composition (113) gives the maximum value Ti for proton and minimum for neutron. 



V. SUMMARY AND CONCLUSION 



With the use of a consistently covariant version of the naive QPM we have studied the spin structure functions 
together with the spin density distributions for the system of quasifree fermions having fixed effective mass xq = m/M 
and the total spin J = 1/2. The main results can be summarized as follows. 

(1) We have shown that the corresponding spin structure functions gi{x) and 52(2;) are mutually connected by a 
simple transformation. At the limit xq — > this transformation is identical to the Wanzura - Wilczek relation for 
the twist-2 term of the g2{x) approximation. At the same time for xo — > the relations for the n — th momenta of 
the structure functions have been obtained. Their particular cases are identical to the known sum rules: Wanzura - 
Wilczek (n = 2,4,6...), Efremov - Leader - Teryaev (n = 1) and Burkhardt - Cottingham (n = 0). Further, we have 
shown how the structure functions are connected with the net spin densities sl(x), st(x) and with the density of the 
total angular momentum j{x). 

(2) The proposed approach has been applied to the description of the proton spin structure with the assumption 
that the proton spin is generated only by the spins and orbital momenta of the valence quarks. Apart from that we 
have assumed that the spin contributions from u and d valence quarks can be defined by the SU(6) symmetry and 
for quark effective mass we used approximation x$ — > 0. We have suggested, how one can in this approach obtain the 
spin functions from the valence quark distributions. Then as an input we have used parameterization of the valence 
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terms resulting from the standard global analysis. On this basis, without any other free parameter, the proton spin 
structure functions and related spin densities have been calculated. Comparison of the obtained structure functions 
gi(x) and gi (x) with the experimental measurement demonstrates that the suggested approach well reproduces the 
basic features of the data on the proton spin structure. 

To conclude, the results presented in this paper and discussion in p2[ suggest, that both the proton structure 
functions g\ and gi have a simple and natural interpretation even in terms of a naive QPM, provided that the model 
is based on a consistently covariant formulation, which takes into account spheric symmetry connected with the 
constraint J = 1/2. This is not satisfied for the standard formulation of QPM, which is based on the one-dimensional 
kinematics related only to the preferred reference system (infinite momentum frame). As a result, there is e.g. the 
known fact, that the function gi{x) has no well-defined meaning in the standard naive QPM. In this case it is just a 
result of the simplified kinematics and not because of an absence of dynamics. 

Acknowledgements I would like to thank Anatoli Efremov and Oleg Teryaev for many useful discussions and 
valuable comments. 



APPENDIX A: PROOF OF THE RELATION (37) 



The integral 



after the substitution 



reads: 



F(x) = / K( Po )6 



Po+Pi 
M 



x I d 3 p; p = \Jm 2 + p 2 



P2 = \l Pq ~ Pi ~ m2 sin f-> Ps = V Po ~ Pi ~ m2 cos f 



F(x) = 2tt / K(po) Po 



+H 



-H 



PO+Pl 

M 



x I dpi 



dp ; H = Jpl - m 2 



For given x and po the inner integral contributes only for 

pi = Mx - po 

in the limits 



ipQ — m 2 < pi < y Po — m 2 . 
One can check, that the last two conditions are compatible only for 

Mx m 2 

Po>Z 



2 2Mx 



It follows, that Eq. (A3) can be simplified 



F(x) = 2irM / K(po)p Q dp 



(Al) 

(A2) 

(A3) 

(A4) 
(A5) 

(A6) 
(A7) 



there are the two roots of x, 



According to the relation (A£) the parameter £ is a function of x with a minimum at xo = m/M, so for the fixed £ 

-, (A8) 



x± 



so the Eq. (A7) can be rewritten: 



M 



M 



2irM / K(po)p Q dpo- 



(A9) 
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Then differentiation in respect to £ gives 



M 



M MVC 2 - m 2 



-2kMK(£)Z 



and this relation can be, with the use of Eq. ( A8), applied to the functions ( 



which with the use of the relation (A6) implies the relation ( p7| 



APPENDIX B: PROOF OF THE RELATIONS (|38|) , (|3g|) 

In the relations ( p2| ) - one can, due to the S— function, make the following substitutions. First, 

Pi = Mx — po 



and then from the relation 



one gets 



Now the relation <\24) can be rewritten: 



2222 
Pt = Po ~ Pi - m 



p 2 T = 2Mxp - (M 2 x 2 + to 2 ) 



gi{x)+g 2 (x) = - 



'Mx Po - (M 2 x 2 + to 2 ) /2\ f po+pi 



xo v-i(x) + 1 H( P0 ) i _ v: : _v "~ i s i - * i d? P 



Po ipo + m) 



M 



In the next step we expand the fraction 



1 1 



E 

3=0 



then the relation (|B4| ) can be expressed in terms of the functions ( p6[ 

1 



gi(x)+g 2 {x) = 



2 i Jl 



xov-i(x) (- x °y e (~ x °y v -^ 



i=2 



2 oo 



(x + x ) V-i{x) - ~-±^-- J2 (-xoY V-i (x) 



Further, from the relation ( |37| ) one obtains 

V-j(x) 



V 2, V T/ , 



\ ,.2.-2 / 



and because 



E 

J=2 



-2xgy V _ 4a; 2 y 2 

2 + ^o/ " (2/ + ; ro) 2 (y 2 + ^) : 



the relation (B6) can be modified 



14 



gi(x)+g 2 (x) = - 



(x + .T ) 



2y T// , . {x + xo) 



2 „1 



4a; y 



o A (y + to) (y 2 + xg) 



V?(i/)dy 



Now, the integration by parts gives 



9i(x)+g 2 (x) = - 



{x + Xo) x* + xJ 



2x (x + xq) 



4t 2 t 2 



2.r 



(t + t ) (x 2 +xl) 



V (x) 



(B9) 



(BIO) 



1 -ix + to) 2 ^-j + ( *+? )2 gffagzfS ^ 

V (2/ 2 + ^) 2 2 ^ (y 2 + ^) 2 (y + T )V 



and one can check, that after some modifications, the corresponding terms ahead of Vq coincide with the functions 
(|l|). The relation @ can be proved by the similar procedure, so we suggest only main steps: 



92{x) 



tVU(t) - V (x) + / H( P0 ) ( Pt-3Mx P0 + (3MV + m i)/2} g ( ^ _ ? 



Po (Po + ™) 



ff2(x) = -- 



o T 2 1 T 2 00 

tV_ x (t) - KoCx) + V (-t ) j V-i(x) + 3- V (-xoY V.j(x) + —V 1 E ^-iO*) 



00 



3=0 



Xq 



3=1 



2x o j=2 



\2 OO 



(2 , r -, ;J r l( ,, f ^^(_ Xo) ^ ]W 



2Xn 



3=2 



52 (x) = - 

Then, the integration by parts gives: 
92(x) = -i 



, / 1 2y 3(t + t ) 

(2t + t ) / ^ — ^ 2 y (y)^ 



2 „i 



Ax 2 Q y 2 



2x1 L (y + T ) 2 (2/ 2 +T 2 ) 



V '(y)dy 



(2.T + T ) 2 

T^ + Xn 



2t 3(t + t ) z 



4t 2 t 2 



2.x; 



(t + t ) (t 2 +Tq) 



1 '(2, + tq) 2 ^ - -$ - 3 ^ 8 ^ ^ Z^L- ) y o(y)dy 
" ' " 2x l (y 2 + xlf(y + x f 



(y 2 + 4) 

and one can check, that the corresponding terms ahead of Vq coincide with the functions (M 



APPENDIX C: PROOF OF THE RELATIONS (|50|) - (|53|) 

First, if we define the functions 



h(y) 



2 2 
V - Xq 

{y + x ) A xo 



h(y) 



y 2 ~Xq 
(y 2 +t^) 2 t ' 



(CI) 



A(x) T 



U(x) = A(x)+a(x) 



2x 2 



a(x) X — To ' '' 1 V ' W T 2 + Tq ' 

then taking into account Eqs. (||3) - (Ell), the Eq. (Ej3) can be rewritten 



(C2) 
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(x + x Q ) 2 fx(y)V (y)dy + (x 2 - x 2 ) f 2 (y)V (y)dy 



(C3) 



-f 3 (x)3(x + x ) 2 f 1 (y)V (y)dy + Mx)(3x 2 + 2xx + x 2 ) f 2 (y)V (y)dy 



ci{x)fi(x)V G (x) + a(x)2 (x + x f / fi(y)V (y)dy - a(x)2 (x 2 + xx + xl) / f2(y)V (y)dy 



-c 2 (x)a(x)V (x) - c 2 (x)3 (x + x ) 2 / fi(y)V (y)dy + c 2 (x) (3x 2 + 2xx Q + x 2 ) j f 2 (y)V (y)dy. 

J X J X 

After differentiating of the l.h.s. one gets 

-2(x + x ) f f 1 (y)V (y)dy + 2x f f 2 (y)V (y)dy 



(C4) 



-3xn 



x 2 — 2xxn — 3x 2 '" 1 



(X - Xq) Jx 



° / fx(y)V (y)dy + 3x 



X 2 — 2xXn — X? r ' 



(x - Xq) Jx 



° I f 2 (y)V (y)dy 



(x + x y h(x) - (x 2 - x 2 ) f 2 (x) + f 3 (x)3 (x + x o y fi(x) - h(x) (3x 2 + 2xx + x 2 ) f 2 {x) V (x) 



2x 2 



= c x (x) 2 V (x) + c 1 (x)2 (x + x Q ) / fi(y)V Q {y)dy - c x (x)2 (x 2 + xx + x 2 a ) j f 2 (y)V (y)dy 

X + Xq 



c 2 {x)2x X 2 X ° 2 Vq(x) - c 2 (x)3 (x + x ) 2 / fi(y)V (y)dy + c 2 (x) (3x 2 + 2xx + xl) / f 2 (y)V (y)dy 



This equality contains the three linearly independent terms 

V (x), f fi(y)V (y)dy, I f 2 (y)V (y)dy 

J x J X 

and comparison of these terms on both the sides gives the equations: 

x + 2xq 



Ci(x)x - C 2 (x) (x - Xq) 



X + X 



(C5) 



(C6) 



2c\{x) — 3c 2 {x) 



2x 2 ~ xxq — 7xq 

(x + X ) (x - Xq) 2 



(C7) 



2ci(x) (x 2 + XXq + Xq) — C2(x) (Sx 2 + 2XXQ + Xq) = 



— x 2 xq — 4xxq 3xq 
(a; - x ) 2 



(C8) 



One can check, that these three equations with the two unknown cx, c 2 are dependent and solvable, giving the solution 
Now, we can solve differential equation (|5fj) for <?i or g 2 . Its homogenous version for gi reads: 

9'i (*) + ( f + tTto\ 9i (x) = 0, (C9) 
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which has the solution 



t + 2t 



Nonhomogeneous term [the part of Eq. (p0|) involving c/2] gives the equation for the function C(x) 



C'{x) 



which has the solution 



, u ^x + 2xq _ x (x 2 + xx a + 4x1) 1 



Xq (x — x ) (x 2 — xf) (x + 2t ) 



92{x) 



X — .Tq 



-92{x) , 



C(x) = ~ x * +2x Q 92 ^ ~ I • / - , -' /) ' /// - 



(CIO) 



(Cll) 



(C12) 



After inserting into Eq. (CIO) one gets the relation (p3|). The inverse relation (p2|) can be proved in a similar way. 



APPENDIX D: PROOF OF THE RELATIONS (65), (fed) 



The relations (B8h and (|3 



92(x) = 



X-Xq 

x x^T4 Voix) 



1 f3 (x + x ) 2 3x 2 + 2xx + t 2 , \ y 2 



(y + xo) (y 2 +x 2 Y 



2.T( 



CjL V (y)dy 



gi(x) +92{x) 



(x + xq) x 2 - x 2 , \y 2 -xf ) 



XXQ 

—5 t> Vq(X) -1/1 , 

x2 + *l h\ (y + xof (y 2 + x 2 f 

can be combined in such a way that the second terms in the integrals cancel: 

(3x 2 + 2xx + x ) gi(x) + 2 (x 2 + xx + xf) 92{x) 



2to 



Vo(y)dy 



(Dl) 



r 1 y 2 _ x 2 

= (x + 2t ) xVo(x) - (x + 2x ) (x + t ) 2 \ -f- ^V (y)dy. 



x (y + x ) 



Further, this equation can be modified 



^-,V (x) - f J^^ Vo (y)dy = ^ + 2 ™° + *° U {x) 



2 (x 2 + xx Q + xf) 

(x + x o y J x (y + t )* (x + 2x Q ) (x + xqY " (x + 2x ) (x + x )' 

which implies a differential equation for Vo(x): 

( x „, A x 2 - xl T _. . / 3x 2 + 2xx + xl . . 2 (t 2 + xx + xf) N 
\(t + t ) J (x + x ) y (x + 2x ) (x + xq) (x + 2x Q ) (x + xq) 

The corresponding homogenous equation 



:92(x) 



gives the solution 

and for C{x) we have the equation 



(x + x ) 

V {x) = C 



V^x) = 



(X + .T ) 



x ,. . / 3t 2 + 2tt +t^ 2 (t 2 + tt + t^) 

2 C ( X > = 7 T^fllW + , ' w -2 92{X) 



(x + 2x ) (x + XqY 



(t + 2t ) (t + T ) 



(D2) 



(D3) 



(D4) 



(D5) 
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The solution reads: 



w s n( x 3a; 2 + 2xx + x 2 2 (x 2 + za; + x 2 ) 

W = = (x + 2x )x gl(a:) + (x + 2x )x g2(x) (D6) 



1 (3y 2 + 2^Q+X 2 ) (y- X °) gi ( y ) dy + ^ 2 + + X ^( y ~ XQ \ g^y)dy 



(y + x ) (y + 2x ) y 2 J x (y + x ) (y + 2x ) y 

After inserting from Eq. (f52"|) one gets 



3x 2 + 2x.t + x 2 , 2 (x 2 + xx + x 2 ,) / a; - x x(x + 2x ) f 1 y 2 - x 2 \ 

MX) = (* + 2x )x glW + (x + 2x )x + l^rX -J^r-ft^J 

/- 1 (3y 2 + 2yx + x 2 ) (y - x ) 
+ L (y + X o)(y + 2x )y 2 ^ 

1 2 (y 2 + yx Q + x 2 ) (y - x ) ( y - x y (y + 2x ) f 1 z 2 - x\ \ 

(y + X0 )(y + 2x )y 2 {-—^ y) + l^ofi —^ z)dz ) dy > 

The double integral is calculated according to the formula 

I" a(y)( [ b(z)dz)dy = [ (A(y) - A(x)) b(y)dy; A'(x)=a(x), (D7) 

x \J y J J x 

then after collecting the corresponding terms wit h g\ we obtain with the use of Eq. ( |64| ) the relation (|65|). A similar 



procedure with inserting from Eq. (p3[) into Eq. (D6) gives the relation (pq 



APPENDIX E: PROOF OF THE RELATIONS @, (frt|) 

In the relation ( |69| ) one can, due to the S function, make the substitution pi = Mx — pq. Then, using the definition 
(j36|), one obtains 

s L (x) = gi (x) + ^V Q (x) - ^xV^(x). (El) 



Further, the relation (37) implies 



2t I" 1 ii 2 — r 2 

v -&) = - 2 / 2 \, Vo(y)dy, (E2) 

x + x o Jx {y z + xf j ) 



which after inserting into Eq. (El) gives 



s L (x) = gl (x) - f X ° V (x) + x f V l X l V Q {y)dy. (E3) 
i\x -\-x ) j x (jj^ + xfi) 

After inserting Vq = 2j(x) and using the relation ( |65| ) we obtain 

, , 1 x Xr, , . x Xr, 3x ~\~ 2xxo 3xn f y Xr. . . 

Si( ^2— ^^(x 2 -^ 2 ) (x + , o) 2 X — ^ 



1)7. ln i^^fJ— 9l(y)dy+a; X (^!)7 gl( ^ 
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1 y 2 -xl 3y 2 + 2yx + 3x 2 f 1 z 2 

(y 2 + xlf (y + x f 



—gi(z)dzdy 



+2x 



1 v 2 



(u 2 



xlf J„ 



In 



y(z + x ) \ £ 
z(y + x a ) 2 I 



—gi(z)dzdy. 



Then calculation of the double integrals with the use of relation (|D7f) and collecting the corresponding terms with gi 
give the relation (|70|). 

Further, one can insert gi from the relation (^3|) into relation ((70|), then a similar procedure with double integrals 
gives the relation (f7l|). 
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